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Within the framework of shallow-water magnetohydrodynamics, we investigate the linear 
instability of horizontal shear flows, influenced by an aligned magnetic field and strat¬ 
ification. Various classical instability results, such as Holland’s growth rate bound and 
Howard’s semi-circle theorem, are extended to this shallow-water system for quite gen¬ 
eral profiles. Two specific piecewise-constant velocity profiles, the vortex sheet and the 
rectangular jet, are studied analytically and asymptotically; it is found that the magnetic 
field and stratification (as measured by the Froude number) are generally both stabilising, 
but weak instabilities can be found at arbitrarily large Froude number. Numerical so¬ 
lutions are computed for corresponding smooth velocity profiles, the hyperbolic-tangent 
shear layer and the Bickley jet, for a uniform background field. A generalisation of the 


long-wave asymptotic analysis of Drazin & Howard (1962) is employed in order to under¬ 


stand the instability characteristics for both profiles. For the shear layer, the mechanism 
underlying the primary instability is interpreted in terms of counter-propagating Rossby 
waves, thereby allowing an explication of the stabilising effects of the magnetic field and 
stratification. 
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1. Introduction 

The interaction of horizontal shear flows and magnetic fields in stably stratified layers 
is central to many problems in astrophysical fluid dynamics — involving, for example, 
planetary interiors, stellar radiative zones and accretion discs. An important example 
of such a flow, which has received considerable attention recently, is that of the solar 


tacliocline (see Hughes, Rosner & Weiss 2007). The tachocline, discovered via helioseis- 


mic observations, is a thin layer in the Sun, extending downwards from the (neutrally 
stable) base of the convective zone to the (stably stratified) top of the radiative interior, 
characterised by radial velocity shear and also planetary scale horizontal shears, associ¬ 
ated with the equator to pole differential rotation of the Sun. Most models of the solar 
dynamo invoke the tachocline as the site for the storage and generation of the Sun’s 
strong, predominantly toroidal magnetic field. 

Here we are interested in the stability of such shear flows, and how this depends upon 
the velocity profile, magnetic field strength, and stratification. Specifically, we consider 
the linear stability of a steady parallel flow and aligned magnetic field, both sheared in 
the horizontal cross-stream direction, in the inviscid and perfectly conducting limit. In 
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this first study, we consider the case where there is no background rotation. The nonlinear 
regime of such instabilities typically leads to turbulent flows; these may be important for 
dynamo action, through some mean-field a-effect, and also for the transport of mass and 
momentum, which can feed back on the large-scale flow. 

It is possible to examine the stability of such flows in a continuously stratified three- 
dimensional setting (e.g., Miura fe Pritchett||1982 Cally 2003). However, here we adopt 
the alternative approach of considering the dynamics of a thin fluid layer under the 
shallow-water approximation, which is valid when the horizontal length scale of the mo¬ 
tion is long compared with the depth of the fluid layer, as is typically the case in large-scale 
astrophysical flows. This leads to a set of two-dimensional partial differential equations, 
with no explicit dependence on the vertical co-ordinate, which offers a considerable math¬ 
ematical simplification. Such shallow-water equations capture the fundamental dynamics 
of density stratification, including gravity waves, and allow the interaction of stratifi¬ 
cation with horizontal shear flows and magnetic fields to be analysed in the simplest 
possible setting. 

It should be noted that hydrodynamic shallow-water models, which date back to 
Laplace, are derived by considering a thin fluid layer of constant density bounded below 
by a rigid medium and above by a fluid of negligible inertia (e.g., Vallis[2006, §3.1). The 
corresponding reduction for electrically conducting fluids — the shallow-water magneto- 


hydrodynamic (SWMHD) equations of Gilman (2000) - additionally requires the fluid 

layer to be a perfect conductor, and to be bounded above and below by perfect conduc¬ 
tors. There are few direct astrophysical analogues for such a configuration. However, we 
can borrow an important idea from planetary atmospheric dynamics, where the hydro- 
dynamic shallow-water equations are widely used to understand waves and instabilities 
in a continuously stratified atmospheric layer. This is justified because there is a for¬ 
mal mathematical analogy between the linearised equations in the two systems, provided 
the layer depth in the shallow-water model is taken to be a so-called equivalent depth 
(e.g., |GflI|[l982l §6.11), so that the shallow-water gravity wave speed (in the horizontal) 
matches that of (the fastest) gravity waves in a continuously stratified layer. We have 
this analogy in mind throughout this study. 

The SWMHD equations have been studied widely in recent years. They have been 
shown to possess a hyperbolic as well as a Hamiltonian structure (De Sterck 2001 


Dellar 


2002), and to support wave motions such as inertia-gravity waves and Alfven waves 
^ Schecter et aL|200l |Zaqarashivili et cd.|2008 Heng fe Spitkovsky|2009 ). As reviewed by 
Gilman & Cally (2007), they have also been used to study the linear instability of shear 


flows in spherical geometry, often with modelling tachocline instabilities in mind. These 
previous studies considered basic states that were functions only of latitude, investigating 
the dependence of the instabilities on the strength and spatial structure of the magnetic 
field and on a (reduced) gravity parameter ( Gilman fe Dikpati|2002 Dikpati et a7|2003 ). 

In contrast to previous investigations of the instabilities of shear flows in SWMHD, 
which focused on global instabilities in spherical geometry, here we concentrate on local 
instabilities, with the aim of examining the linear instability problem in a wider con¬ 
text; for this, we consider the problem in planar geometry. We first derive some growth 
rate bounds and stability criteria, valid for general basic states. We then study how the 
instability characteristics of prototypical flows are modified by the combined action of 
magnetic fields and stratification, which, in isolation, are generally thought to be sta¬ 
bilising. The corresponding hydrodynamic problem has a long history, dating back to 
Rayleigh, and we are able to draw upon ideas and methods from a substantial literature 
(e.g., Drazin fe Reid|[l98l Vallis|[2056 ). 

We start, in §2, by formulating the linear instability problem for plane-parallel basic 
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states with the flow and field dependent on the cross-stream direction. In §3, we de¬ 
rive extensions of classical growth rate bounds, semi-circle theorems, stability criteria, 
and parity results for modal solutions. In §4, instabilities of piecewise-constant velocity 
profiles (the vortex sheet and the rectangular jet) with a uniform magnetic field are ex¬ 
amined analytically. Analogous smooth profiles (the hyperbolic-tangent shear layer and 
the Bickley jet) are studied in §5, both numerically and asymptotically, via a general¬ 


isation of the long-wave analysis of Drazin & Howard (1962). The primary instability 


mechanism for shear layers is interpreted in terms of counter-propagating Rossby waves. 
The results are discussed in §6. 


2. Mathematical formulation 

2.1. Governing equations 

We consider a thin layer of perfectly electrically conducting fluid moving under the 
influence of gravity. We use a Cartesian geometry, with horizontal coordinates x and y, 
and an upwards pointing coordinate z. At time f, the fluid, which is taken to be inviscid 
and of constant density p, has a free surface at z = h{x , y, t ) and is bounded below by a 
rigid impermeable boundary at z = —H(x,y). 

We consider motions with a characteristic horizontal length scale Lq that is long com¬ 
pared with a characteristic layer depth H 0 . One can then make a shallow-water reduction 
in which the vertical momentum balance is taken to be magneto-hydrostatic, and for 
which the horizontal velocity u and horizontal magnetic field B are independent of z. 
When the bottom boundary is perfectly electrically conducting (and is thus a magnetic 
field line) and the free surface remains a field line, the magnetic shallow-water equations 
of Gilman (2000) are obtained. These are an extension of the classical shallow-water 


equations of geophysical fluid dynamics. 

We use these equations in non-dimensional form. We denote the characteristic hori¬ 
zontal velocity of the basic state by C/o, and the characteristic magnetic field strength 
by B 0 . Non-dimensionalising x and y by Lq, t by the advective time-scale Lq/Uq, H by 
Ho , h by Uo/g (where g is the acceleration due to gravity), velocity by Uo , and magnetic 
field by B 0 , the SWMHD equations are 

Vit = —V/i + M 2 B ■ X7B, (2.1 a) 


dt 

dB 

dt 


+ u ■ 


u • VB = B • Vit, 


, dh 


F 2 Q t +V-(( H + F 2 h)u)=0, 


(2-1 b) 
(2.1c) 


where F = U 0 /\/gH 0 and M = {B 0 /y/JIp)/U 0 , with p being the permeability of the fluid. 
In addition to (2.1 z-c), the shallow-water reduction implies 


V • ((H + F 2 h)B ) = 0. 


( 2 . 2 ) 


However, (2.2) need not be considered explicitly; if it is satisfied at some initial time, 


then (2.1 z-c) guarantee that it remains satisfied for all time. 

The system has two non-dimensional parameters. The Froude number F is the ratio of 
the characteristic horizontal velocity of the basic state to the gravity wave speed \fgHo 
(and is related to the reduced gravity parameter G of Gilman & Dikpati (J2002| via 
G = F ~ 2 ). The parameter M is the ratio of the Alfven wave speed Bo/^/JIp to the 
characteristic horizontal velocity of the basic state. When H is constant and F —> 0, 
(12.1 cl) and 12.21) become V • u = 0 and V • B = 0 respectively, and we recover the 
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equations for two-dimensional incompressible magnetohydrodynamics, with h playing 
the role of pressure. When M — )■ 0, (2.16) decouples from (2.1r-c), and we recover the 
hydrodynamic shallow-water equations; these have a well-known correspondence with 
two-dimensional compressible hydrodynamics (e.g., Vallis 2006 §3.1), which we exploit 
from time to time. 

As an example of astrophysical parameter values, we estimate F and M in the solar 
tacliocline, using data from Gough (2007). We set Uq to be the equator to pole difference 
in the zonal velocity, implying /7o ~ 500ms -1 . There is considerable uncertainty in the 
strength of the magnetic field in the tachocline (Hughes et al. 2007), although a likely 
range is 10 3 G < Bq < 10 5 G. Then, taking p = 210kgm -3 , we find 0.01 < M < 1. 
To estimate F, we must choose a gravity wave speed \J gH o for the layer. One means of 
doing this is to take H 0 to be the depth of the tachocline and to interpret g as a reduced 
gravity, accounting for the fractional density difference of the overlying fluid, as in|DikpatI| 


& Gilman (2001). However, here we pursue the analogy between shallow-water flows and 
those of a continuously stratified layer with buoyancy frequency N and depth Hi, and 
ch oose V g Hp to be the speed of the fastest gravity wave in such a layer, which is NHi/tt 
(Gill 1982 §6.11). Taking Hi ss 2 x 10 7 m (i.e. O.O3i? 0 , where Rq is the solar radius) and 
N S3 8 x 10 -4 s -1 , which are bulk values that might describe a mode spanning the entire 
tachocline, gives a gravity wave speed NHi/tt = y/gH 0 ss 5000ms -1 , corresponding to 
an equivalent depth Hq ~ 50 km (taking g ss 540 ms -2 ). Again taking Uq ss 500 ms -1 , 
we thus estimate F ss 0.1, although it is clear that F would be somewhat smaller or 
larger if one considered motions towards the top of the radiative zone (with stronger 
stratification) or towards the base of the convection zone (with weaker stratification). 


2.2. The linear instability problem 

Above a topography of the form H = H(y), we consider a basic state h = 0, u = U(y)e x 
and B = B(y)e x , so that the magnetic field is initially aligned with the flow. We then 
consider perturbations in h, u = ( u , v) and B = ( b x , b y ) to this state of the form 

£(x, y, t ) = R e{i(y) exp (ia{x - ct ))}, (2.3) 


where a is the (real) wavenumber and c is the (complex) phase speed. Dropping the 
hatted notation, the linear evolution is described by 


d 


d 


dt 


dx 

d 


dh 


— + U— u + U'v = — — + M 2 B—^ + B'b. 


/- + U-/- ) v = + M’B^, 


dt 

d 

dx 

d 


dx 


dx 

dh 


. db a 


dx 


v I 5 


dy 

du 


dx 


ft + U TH ) b * + B ' v ~ B HZ + U ' b v ’ 


dt 


dx 


F 2 — + U— )h + H[— + —)+H , v = 0, 


dt + U dx 


du dv 


dx 

by ~ dx' 


dx dy 


(2.4a) 
(2.46) 
(2.4c) 
(2.4 d) 
(2.4e) 


where a prime denotes differentiation. Eliminating for v, we obtain 

/ S 2 {Hv)' V / a 2 S 2 U' ( S 2 y Q'S 2 \ 

\H(U — c) 2 K 2 J [H(U-c) 2 H(U -c)\(U- c) 2 K 2 ) + {U - cfK 2 ) 


Hv = 0, 
(2.5) 
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where Q = —U'/H is the background potential vorticity, and 

S 2 (y) = ( U(y ) - c) 2 - M 2 B 2 (y ), I< 2 (y) = 1 - F 2 S 2 (y). 


( 2 . 6 ) 


Following Howard (1961), under the transformation Hv = (U — c)G, equation (2.51 
becomes 

S 2 G' V a 2 S ? 

K 2 H ) H 

We shall use this more compact form for the remainder of this study. In the non-magnetic 


G = 0. 


(2.7) 


shallow-water limit (M = 0), (2.5) reduces to equation (3.4) of Balmforth (1999). In the 
two-dimensional incompressible magnetohydrodynamic limit (F = 0 and Ft = 1), (2.7) 
reduces to equation (3.5) of Hughes & Tobias (2001). 

We shall consider (2.7) in either an unbounded domain, for which |G| —>■ 0 as \y\ —► oo, 


or in a bounded domain with rigid side walls, where G = 0 and hence b y = 0 via 


(2Ad). Either way, for given real a, (2.7) is then an eigenvalue problem for the unknown 
phase speed c = c r + ic*. We will focus on instabilities, i.e. c,; ^ 0, in which case (2.7) 
has no singularities for real values of y. Since the transformation a —> —a leaves (2.7) 


unchanged, we may take a ^ 0 without loss of generality. Instability then occurs if c; > 0, 
with growth rate aCi. 


3. General theorems 

In this section we derive three results that hold for general shear flows U(y): two provide 
bounds on the growth rate of any instability, whereas the third concerns implications of 
the parity of the basic state flow. 


3.1. Growth rate bound 

A bound on the instability growth rate may be obtained by calculating the rate of change 
of the total disturbance energy using the combination 


Hu* x (2.4a) + Hu* x (2.46) + {M z Hb* x ) x (2.4c) + (M 2 Hb* y ) x (2Ad) + h x (2.4e), 


where * denotes complex conjugate. On adopting the form ( |2.3[ ) for the perturbations, 
the real part of this expression gives (on dropping hats) 


aa (H (M 2 + \v\ 2 + M 2 |6,| 2 + M 2 \b v \ 2 ) + F 2 \h \ 2 ) = 

- Re (HU 1 (vu* - M 2 b* x b y ) + M 2 HB' (vb* - u*b y )) - Re|(M*). 


On integrating over the y domain, employing the boundary condition on v , and manip¬ 
ulating the remaining terms on the right hand side using ±2Re(pg*) < |p| 2 + \q \ 2 , we 
obtain the following bound on the growth rate: 

ac i ^ -(max \U'\ + Mmax \B'\). (3.2) 


In the absence of magnetic field, this reduces to the well-known bound in hydrodynamics 


( H0iland||1953 Howard ||~1 961 1 . 


3.2. Semi-circle theorems 


In a classic paper, Howard (1961) proved that for incompressible hydrodynamic parallel 
shear flows, the wave speed c of any unstable mode must lie within a semi-circle in the 
complex plane determined by properties of the basic state flow. Subsequently, semi-circle 
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theorems have been derived for several other hydrodynamical and hydromagnetic systems 


(e-g-, 

Codings & Grimshaw 1980 

Hayashi & Young||1987 

Shivamoggi & Debnath 1987 

{Hughes & Tobias 

2001 

). In a similar manner, a semi-circle theorem may be derived for 


the SWMHD system. 


Multiplying equation (2.7) by G*, integrating over y and using the boundary condition 

S 2 \G\ 2 


on v (and hence G) gives the relation 

f S 2 | G' 12 


J K 2 H 

The imaginary part of (|3.3[) gives 


d y + a 


■ J {U - c r )x d y = 0, where 


H 


X = 


dy = 0. 


| G' 


H\K\ 


(3.3) 


4 ,a 2 ^-^0. (3.4) 


Equation (3.4) immediately yields Rayleigh’s result that for unstable modes (c, > 0), c r 
lies in the range of U (i.e. C/ m i n ^ c r ^ t/ max , where the subscripts ‘min’ and ‘max’ refer 
to the minimum and maximum values across the domain). 

On using equation (3.4), the real part of (3.3) gives 

(c 2 + c 2 ) J X dy = j X (U 2 -M 2 B 2 ) dy - F 2 J | G" | 2 dz/, (3.5) 

which implies that 

0 ^ (c 2 + c 2 ) J x dy B (U 2 - M 2 B 2 )_J x dy. (3.6) 

This gives the first semi-circle bound: the complex wave speed c of an unstable eigen¬ 
function must lie within the region defined by 

c 2 + c 2 < (U 2 -M 2 B 2 ) max . (3.7) 

The second semi-circle bound is obtained, in the standard manner, from the inequality 


0 ^ f(U — t / max )(U — Brain) X dy ■ Substituting from (3.4) and deriving an inequality from 
(|3.5|) leads to the expression 


0 JS (c 2 + C 2 — (Um in + Umax)c r + GninKnax + M 2 (B 2 ) m in) J X 


dy, 


(3.8) 


which gives the second semi-circle bound: the speed c of an unstable eigenfunction must 
lie within the region defined by 


Cr* 


U m in + Ur r 


+ c i ^ 


B i 11 ;i.x Ur r 


— M 2 {B 2 ) T 


(3.9) 


Thus, taking these results together, the eigenvalue c of an unstable mode must lie within 


the intersection of the two semi-circles defined by (3.7) and (3.9). In the absence of 


magnetic field, semi-circle (3.9) lies wholly within semi-circle (3.7), and we recover the 


well-known result of Howard (19611. However, as observed by Hughes & Tobias (2001), 


who considered the stability of aligned fields and flows in incompressible MHD, for non¬ 
zero magnetic field there is the possibility of the two semi-circles overlapping, being 
disjoint, or indeed ceasing to exist; thus, in addition to giving eigenvalue bounds for 


unstable modes, these results also provide sufficient conditions for stability. From (3.7) 


and (|3.9|) it therefore follows that the basic state is linearly stable if any one of the 

(3.10) 


following three conditions is satisfied: 

M\B\ ^ | U | everywhere in the domain; 
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M\B\r 




| Gmax - U n 


(3.11) 


U„ 


U r r 



~U„ 


+ M 2 (S 2 ) min ^ J(U 2 -M 2 B 2 ) 


(3.12) 


These results are equivalent to those given by Hughes & Tobias (2001) for incompressible 

MHD. 

A drawback of the above approach is that the bounds do not contain the Froude 
number F. Although it is possible to introduce F into the semi-circle bounds using 


similar manipulations to that employed by Pedlosky (1964), as shown by Mak (2013) 


this does not sharpen the bound and we thus omit it. 


3.3. Consequences of basic state parity 

For the hydrodynamic case, it can be shown that symmetries of the basic state lead to 
symmetries in the stability problem (H oward|[l963 ). These results may be generalised to 
SWMHD if we make the further assumptions that B 2 (y) and H(y) are even functions 
about y = 0. 


We first consider the case when U(y) is odd about y = 0. Equation (2.7) is unchanged 
under c —► — c and G(y) —> G{—y). Since the equation is also unchanged under c —> c* 
and G —> G*, it follows that an eigenfunction with eigenvalue c = c r + Ci must be 
accompanied by eigenfunctions with c = ±c r ± icj. Thus unstable solutions either have 
c r = 0 or are a pair of counter-propagating waves with the same phase speed. As argued 


by Howard (1963), the symmetry in the basic state implies that there is no preferred 


direction for wave propagation, consistent with the form of the eigenvalues. 
Now consider the case when U(y) is even about y = 0. Then 


°e(y) = \ ( G{y ) + G{-y)) and 


Go(y) = \ ( G(y) G{-y )) 


(3.13) 


are also eigenfunctions of (2.7). Following Drazin & Howard (1966), if we now take G 0 


multiplied by (2.7) with G = G e and subtract this from G e multiplied by (2.7) with 
G = G 0 , integrating over —L y ^ y ^ L y gives 


W(G e , G 0 ) = [G' e G 0 - G’ 0 G e ] + _ L L y y = constant = 0, (3.14) 

owing to the imposed boundary conditions on the eigenfunction. The vanishing of the 
Wronskian W implies that the functions G e and G 0 are linearly dependent throughout 
the domain, which is possible only if one of them is identically zero. Thus an unstable 
eigenfunction corresponding to a particular eigenvalue is either an even or odd function 
about y = 0. 


4. Piecewise-constant profiles: vortex sheet and rectangular jet 

We now consider some simple flow configurations for which the eigenvalue problem (2.71 
can be reduced to an algebraic equation for c, from which the conditions for stability can 
be readily determined. To do this, we take H = 1 (no topography) and B = 1 (a uniform 
magnetic field). We seek solutions of (2.7) in an unbounded domain, with 

|G| —> 0 as \y\ —> oo. 


(4.1) 


We consider velocity profiles U(y) that are piecewise constant. If U(y) is discontinuous 
at y = yo, then the eigenfunction G must satisfy two jump conditions at y = yg. In the 
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usual way, the (linearised) kinematic boundary condition implies 


U-. 


Vo 


Vo 


= \G]f = 0. 

VO 


(4.2 a) 


The pressure (or free surface displacement) is also continuous at y = yo. The correspond¬ 
ing condition on G is most easily derived by integrating (2.7) across y = yo, yielding 


s 2 r , 


vi 


= 0. 


(4.26) 


-‘Vo 

4.1. Vortex sheet 

We first consider the velocity profile 


U{y) = 


+1) 


y > o, 
y < o. 


Then, for y ^ 0, (2.7) becomes G" — a 2 K 2 G = 0. Using (4.1) and (4.2a), we thus find 


G(y) = 


exp (-otK+y ), y > 0, 
exp (+aK-y), y< 0, 


where 


K± = i/l — U 2 ((1 =F c) 2 — M 2 ), Re (A±) > 0. 


(4.3) 


(4.4) 


The second jump condition (4.26) then implies an eigenvalue relation for c: 


= 0 . 


(4.5) 


(1 — c) 2 — M 2 (1 + c) 2 — M 2 
K + + K- 

Note that c is independent of the wavenumber a, so any unstable mode with a > 0 has 
an unbounded growth rate as a —> oo. This is an artefact of considering ideal fluids; 
viscosity will preferentially suppress small scales and remove this unphysical behaviour. 

There are several special cases. When F = M = 0, we recover the classical Kelvin- 
Helmholtz instability with c = ±i. When F = 0 but M ^ 0, ( |4.5| ) reduces to the 
incompressible MHD case of Michael (19551, with c 2 = —(1 — A/ 2 ); thus, the Kelvin 


Helmholtz instability is stabilised when M ^ 1, since the disturbance has to do work 


to bend the field lines. When M = 0 but F ^ 0, (4.5) gives the classical hydrodynamic 


shallow-water dispersion relation, which is analogous to that of two-dimensional com¬ 
pressible hydrodynamics. The Kelvin Helmholtz instability is stabilised when F V y/2 
(Miles][l958 Bazdenkov & Pogutse|[l983 ), since the disturbance has to do work to move 
the free surface against gravity. Thus, increasing F or M in the absence of the other is 
stabilising. 


In the general case where F and M are both non-zero, (4.5) can be rearranged and 
squared to yield a quartic equation for c: 


U 2 c 4 - 2 (1 + F 2 (M 2 + 1)) c 2 + ( M 2 - 1) (2 + F 2 ( M 2 - l)) = 0. 


(4.6) 


Here we have ignored the degenerate case with c = 0, which is a solution of (14. 51) when 


M = 1. By comparing solutions of (4.61 with those of (4.5) found using a Newton iteration 


method, we find that only two roots of (4.6) also satisfy (4.5): these are c = ±c v , where 

( Vl + 4F 2 + 4F 4 M 2 - (1 + F 2 + F 2 M 2 ) \ , N 

Cv = 1 1 -^-- ■ (4.7) 


F 2 
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A contour plot of Im(c v ) is shown in figure [I] From (4.7), there is instability only if 

M < 1 


and 


F < 


1 - M 2 ' 


(4.8) 


Although increasing F is always stabilising at fixed M. the critical value of F above which 
the flow is stable increases as M increases towards 1. Thus, although magnetic field and 
free-surface effects are stabilising in isolation, together they can lead to instabilities at 
arbitrarily large values of F, provided 


1- — <M 2 <1. 
t z 


(4.9) 


Using an asymptotic analysis, it is possible to investigate these instabilities further at 


large F and with M just smaller than unity. Rewriting (4.7) in terms of 1 — M 2 and 
expanding for |1 — M 2 | <C 1, we obtain 


1 — M 2 2F 6 (1 - M 2 ) 2 \ 1/2 
1 + 2 F 2 (1 + 2F 2 ) 3 


11 — M 2 1 < 1, 


(4.10) 


where terms of O ((1 — M 2 ) 3 ) have been neglected. When F = 0(1), the first term on 
the right-hand side of ( ]4.10[ ) dominates. However, in the regime of interest (4.9) with 
F 2 ~ (1 — AT 2 ) -1 1, the two terms on the right-hand side of (4.10) have the same 

order of magnitude, and instead we obtain 


1 - M 2 (1 - Af 2 ) 2 


2 F 2 


1/2 


as F~ 2 ~ (1 - A/ 2 ) -+ 0. 


(4.11) 


This simple formula is consistent with both stability boundaries in (4.8 1 , and, as shown in 


figure [2} clo sely predicts c, in this weak instability regime, even when F is of order unity. 
Using (4.11), it is straightforward to show that Im(c v ) is maximised when A/ 2 = 1 — 1/F 2 , 
with c v ~ i/(2F 2 ), so that the growth rate of the most unstable mode decays like F~ 2 
in this regime. 
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Figure 2. The weak insta bility regime of the vortex sheet for (a) F = 2, (b) F = 5, as 
determined directly from (4.TI (crosses) and from the asymptotic result (4.111 (line). 


4.2. Rectangular jet 

We now consider the top-hat velocity profile 


U(y) = 


i, \y\ < i, 

0, \y\ > 1. 


Then, (2.71 and (4.1) imply 


A + exp(-aK 0 (y - 1)), y > +1, 

G = { A e cosh(aK^y) + A a sinh(a/vij/), \y\ < 1, 

A- exp (+aK 0 (y + 1)), y < -1, 

for some A + , A-, A e and A a , where 

K 0 = \/l — F 2 (c 2 — M 2 ), with Re(K 0 ) > 0 for bounded solutions, 

K\ = \J 1 - F 2 ((1 - c) 2 - M 2 ), with - ^ < arg(ATi) < 

Here we follow Rayleigh’s formulation (Drazin & Reid 1981| ) and consider eigenfunctions 
that are either even or odd. For the even mode, we set A 0 = 0, A + = A_ and write 
c = c e . Then (4.2 ab) and (4.13) give 


(4.12) 


(4.13) 


(4.14a) 

(4.146) 


- M 2 (1 - c e ) 2 - M 2 


+ 


tanh(aA'i) = 0. 


(4.15) 


K 0 K 1 

For the odd mode, we set A e = 0, A + = —A _ and write c= c 0 . Then ( 4.2a| 6) and (4.131 
give 

c 2 - M 2 (1 - c 0 ) 2 - M 2 


+ 


coth(aAd) = 0. 


(4.16) 


A„ K\ 

In contrast to the vortex sheet dispersion relation ( |4.5[ ), here c depends upon a. 

Two special cases may be solved analytically. When F = M = 0, so that Ao,i = 1, 
expressions (4.15) and (4.16) yield unstable modes with 


T + is/r 
l + T ’ 


i + iVr 
1 +t ’ 


where T = tanlra, 


(4.17) 


so that the flow is unstable for all a, with Im(c) approaching a maximum value of 1/2 
as a —> oo (Rayleigh 18781. When F = 0 but M 0, so that K 01 = 1 again, (4.15) and 
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(4.16) yield 


c P = cl°) = 


T + \\JT — M 2 (l + T) 2 
1 +T 


c = c (0) = 

'-'O '-'n 


1 + i\jT — M 2 (l + T) 2 
1 + T 


(4.18) 


(Gedzelman (1973) gave formulae for these two cases (his (3.4) and (3.5)), but both are 
missing factors of 2.) Both modes are stable for all a when M ^ 1/2; otherwise, both 
modes are unstable (with the same growth rate) provided that 


VT l - 2 M 2 -Vl- 4 M 2 

M < Mi 0 ) = --— •£=> a > ai°) = tanh 1 - 


1 + T 


2 M 2 


(4.19) 


Thus, the magnetic field introduces a long-wave cutoff. 

The other special case with M = 0 but F / 0 cannot be solved analytically. Then, 


(4.15) and (4.16) are equivalent to expressions given by Gill (1965), who considered cor¬ 


responding instabilities for two-dimensional compressible hydrodynamics. An important 
property is that, at fixed a and F, there can be a large number of unstable modes for 


sufficiently large F, which are often interpreted in terms of over-reflection (Takehiro & 


Hayashi 1992). This is in contrast to the behaviour at F = 0, where just two modes exist 


(one even and one odd), as described by (4.17). 


We generalise these results to cases with F ^ 0 and M ^ 0 by solving (4.15) and 


(4.16) numerically using a Newton iteration. We start by limiting attention to the smooth 
extensions of the unstable even and odd modes (4.17) that exist at F = 0 and M = 0, 


and then tracking these modes over (M, F) space; the existence of other unstable modes 
at high F is discussed in Q4.2.5 Contours of Im(c) obtained using this tracking approach 
are shown in figure |3j for a wide range of F, M and a. We note that (i) neither mode 
is unstable for M > 1/2; (ii) for small a, the even mode is stabilised at values of M 
considerably less than 1/2; (iii) for small a and small F, the odd mode is stabilised at 
values of M considerably less than 1/2; (iv) both modes remain unstable for large F\ 
(v) when both modes are unstable, the even mode is generally more unstable than the 
odd mode; (vi) at large a , the even and odd modes lead to instabilities of comparable 
strength, which mimic that of the vortex sheet (cf. figure [l]). We now use asymptotic 
analyses to describe this behaviour in more detail. 

4.2.1. Instability at large a 

When M and F are of order unity, tanh(aA'i) ~ 1 when a>l (provided aRe(A'i) 3> 
1), so that both (4.15) and ( |4.16 l may be approximated by 

(c + 1/2) 2 — M 2 (c — 1/2) 2 — M 2 


= 0 , 


(4.20) 


■>Jl - F 2 ((c+ 1/2) 2 - M 2 ) ■ S /1-F 2 ((d-1/2) 2 ~M 2 ) 

where c = c e — 1/2 or c = c Q — 1/2. This dispersion relation is similar in form to that for 
the vortex sheet ( |4.5| ), and a solution may be found in the same way: 

/ Vl + F 2 + F 4 M 2 - (1 + F 2 /4 + F 2 M 2 ) \ 

C = 1 I- Y 2 -) ' (4 ' 21) 

Physically, a sufficiently localised short-wave disturbance sees only one flank of the jet, 


and thus a vortex sheet instability is obtained to a first approximation. From (4.21), 
there is instability only when 


M < - and 
2 


F < 


1 - 4 M 2 ‘ 


(4.22) 
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Im(c) 



0 0.1 0.2 0.3 0.4 0.5 


E = O = 0 0> E >0 


E > O > 0 O > E = 0 


even 


odd 



oo 

o 


a 


cs 

o 


a 


o 

o 
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Figure 3. Instability o f the rectangular jet. Smooth extens ions a t fixe d a fr om ( M,F ) = (0,0) 
of the unstable modes (4.171, computed numerically using (4.151 and (4.161. Note the different 
choice of F-axis used in the top panels. Left: contours of Im(c) for the even mode. Centre: 
contours of Im(c) for the odd mode. The black contours are in intervals of 0.025 from 0.025 to 
0.5; weak instabilities are shown by the grey contours at 0, 0.001, 0.005 and 0.01. The a 1 
stability bo undar ies ( 4.22|) are shown as dashed lines in panels (a,b). The a -C 1 stability 
boundaries (4.271 and (14.30 1 are shown as dashed lines in panels (g. j) and (h, k) respectively. 
Right: regime diagrams deduced from the left and centre panels, comparing the growth rate of 
the even mode ( E ) with that of the odd mode (0). No regions were found with E > O = 0. 
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As shown by the dashed lines in figures |3ji, b, the conditions (4.22) approximately bound 
the region of strong instabilities in (M, F) space, even when a = 10. However, it is also 
clear that there are additional weak instabilities (shown as grey co ntours ) even when 
F > \/8/(l — 4M 2 ); the nature of these modes will be discussed in ^4.2.5 In contrast, 
the M = 1/2 cutoff is robust. Indeed, as is evident from figures [3|l, e, the magnetic field 
cutoff stays close to M = 1/2 (for all F) even when a « 1. For example, at a = 0.8 and 
F = 0, the cutoff M c « 0.490, using (14.191). 


4.2.2. Instability at small a 

We consider first the even mode. Assuming that K\ remains bounded (which may be 
confirmed a posteriori), tanh(aA'i) ~ a A/ when a < 1, so that (4.15) becomes 

c? - M 2 


y/1- F 2 {c 2 e -M 2 ) 


+ a ((1 — c e ) 2 — M 2 ) =0. 


(4.23) 


Suppose that F 2 = 0(1). If M = 0(1), then c e = ±M at leading order, and the 
next correction in a is also real: no instabilities are predicted, consistent with figure [3]/. 
However, if M 2 ~ a <C 1, then at leading order, (4.23) gives 


±i\/a — M 2 


as 


M 2 


a 


0, A 2 = 0(1), 


(4.24) 


which is consistent with the small M and a limits of (4.18). Thus a weak magnetic field 


reduces the strength of the hydrodynamic instability and eventually suppresses it, with 
instability when 


M < M c = a 1 / 2 , or equivalently when a > a c = M 2 . 


(4.25) 


A dependence on F appears only at the next order in a, which is why the cutoff in 
figure [3]/ is approximately independent of F. 

To capture the destabilising influence of F apparent in figure [ 3 J 7 , the analysis may be 


extended to larger values of F with F 2 ~ a 1 1. In this case, the square root in (|4.231 

also enters the balance at leading order, and we obtain 

1/2 


a 2 F 2 


- M 2 + 


a 4 A 4 


as 


M 2 


F~ 


(4.26) 


which formally reduces to (4.24| when A is of order unity. Again there is a long-wave 


cutoff due to the magnetic field, with instability only when 
a 2 A 2 + \f 4a 2 + a 4 A 4 


M < M c = 


a > a c = 


M 2 


Vl + M 2 F 2 ' 


(4.27) 


This cutoff, which is shown as a dashed contour in figures [3]7, j, captures the destabilising 
influence of A. Indeed, in contrast to the result ( |4.22[ ) for a 1, equation (4.26) does 
not predict stabilisation at large A when a <C 1, consistent with the results shown in 
figures [3j/, j. The accuracy of this dispersion relation and cutoff when a< 1 can be seen 
figure [4j 


In the same a -C 1 limit, the dispersion relation (|4.16|) for the odd mode reduces to 

(4.28) 


a ( c o - M 2 ) 


(1 - c 0 f - M 2 


sj 1 - F 2 (c 2 - M 2 ) 1 - F 2 ((1 - Co ) 2 - M 2 ) 


= 0 . 


Suppose that A 2 = 0(1). If M = 0(1), then c Q is real at the first two orders in a: no 
instabilities are predicted, consistent with figure [3^. However, if M 2 ~ a 1, then at 
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Figure 4. Instability of the rectangular jet at M = 0.1, for (a) F = 0.8 and (b) F = 5. Shown 
are the growth rates of the even and odd modes, computed nu merica lly from (4.151 and ( 4.16 1, 
and the asymptotic predictions w hen a <C 1 for the even mode (4.261 and the odd mode (?!!55). 
In (a), the asymptotic prediction (4.211 for both modes when a 3> 1 is also plotted. 


leading order, (4.28) gives 


1/2 

1 + i [ . = - M z ) as M A 


a 


0, F 2 = 0(1), 1 -F 2 = 0(1). (4.29) 


When F < 1, again there is a cutoff due to the magnetic field, with instability only when 


M < M c = 


a 


1/2 


(1 F 2 ) 1 / 4 ’ or equivalently when a > a c = M 2 \J\ — F 2 . (4.30) 

This cutoff, which is shown as a dashed contour in figures 1 h k, captures the sharp 

1 , which is particularly evident in figure H- However, 
expression (|4.29|) shows that c Q always has a positive imaginary part, so 


destabilising transition as F 
when F > 1 


that there is no cutoff at small M (when a <C 1). This is clear in figures [3^i, k , where 
there are weak instabilities for large F with M close to 1/2. The absence of a cutoff when 
F > 1 is also evident in figure [4f?. 

When F is close to 1, the asymptotics leading to ( |4.29 ) break down, and one must 
instead seek solutions with c = 1 + 0(ct 2 / 5 ). The resulting (quintic) equation does not 
admit solutions in closed form (Mak[2013]), although the special case with F = 1 may 

1 _i_ I 


be solved exactly, yielding c 0 = 1 + ar/ b e 37n / 5 /2 1 / 5 . Thus c 0 is independent of M in this 
regime, consistent with the behaviour shown in the bottom left corners of figures [3[i, k. 

4.2.3. Preferred mode of instability: even versus odd modes 

From figure [3j it can be seen that the even mode is more unstable than the odd mode 
for some parameters but not for others. When a 1, both even and odd modes satisfy 
the same dispersion relation (4.20) to leading order, as seen in figures [3jr, 6; as shown in 
figure [3]::, the even mode is generally more unstable than the odd mode, although this is a 
weak effect. In contrast, when a<Cl, the odd mode is generally more unstable, as shown 
in figure |3j. In particular, when F < 1, the odd mode is more unstable (from (4.24) and 
( |4.29[ )) and is unstable over a larger region of parameter space (from ( |4.25 ) and (4.30)). 
When F > 1, again the odd mode is unstable over a larger region of ( M,F ) space; the 
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M, from (4.30). 


even mode is stabilised according to (4.27), whilst the odd mode has no cutoff at small 


The nature of the transition when a is of order unity can be quantified by performing 
a small F analysis. We thus write c = -t- ^M 1 ) 

for c e and (4.16) for c, 


0{F 4 ) and substitute into (4.15) 


The leading order terms ci°' ) and ci° ■* are given by (|4.18|). In the 


case where there is instability at leading order, i.e. when (4.19) is satisfied, we then find 


Im = a e g, Im = a 0 g, where g(a,M) = 

T + T 2 — a(l — T 2 ) 


T 2 - 6T + 1 
(T+l) 2 
T + T 2 + aT(l-T 2 ) 


-AM 2 


(4.31) 


4(i + t) 2 v / t^Tm^(TTtF’ ° 4(i + t) 2 v / t^Tm^(ITtF’ 


and T = tanha. It is easy to see that a 0 > a e > 0 for all a, so the behaviour is 
determined by the sign of g. More precisely, the even mode is more unstable (at small 
F ) when AM 2 < (6 T — 1 — T 2 )/{T + l) 2 . This is impossible when T < 3 — 2\[2, so the 
odd mode is more unstable at small F when a < 0.1733, for all M. If T > 3 — 2-\/2, then 
the even mode is more unstable at small F when 


,, \J6T — 1 — T 2 , , / 3 - AM 2 - 2V2 - AM 2 

M < M c = ——--—- a < tanh — 


2(1 + T) 


1 + AM 2 


(4.32) 


For example, when a = 0.2, the most unstable mode at small F is even when M < 0.159 
and odd when 0.159 < M < 0.371 (cf. figure |3j); when a = 0.8, the most unstable mode 
is even when M < 0.479 and odd when 0.479 < M < 0.490 (cf. figure [3]/). As a —» oo, 
T —> 1 and M c —»• 1/2, so there is only a vanishingly thin region close to M =1/2 where 
the odd mode is more unstable at small F. 


4.2.4. Maximum growth rate 

At fixed F and M, it is natural to try to find the most unstable wavenumber by 
maximising aci with respect to a. However, this is not guaranteed to lead to a finite 
value of a: for example, when F = 0, (4.18) shows that Ci —> y/l — AM 2 /2 asa-> oo, so 
provided M < 1/2 the growth rate aci increases without bound as a —» oo. The same 
phenomenon persists for small F, as illustrated in figure [l]r at M = 0.1 and F = 0 .8. 
However, when F > ^8/(1 — AM 2 ) and M < A/2, the asymptotic theory of { 4.2.1 for 


a 1 predicts that the flow is stable. In practice, we have seen in figures |3|2, b t hat 
there are weak instabilities for these parameters, but the growth rates decay sufficiently 
rapidly with a that a most unstable mode typically occurs at finite a , as illustrated in 
figure [4f> at M = 0.1 and F = 5. However, we have been unable to find an analytical 
expression for the most unstable wavenumber in this regime. 


4.2.5. Multiple modes of instability at large F 

In addition to the smooth extensions of the unstable even mode and odd mode that 
exist at ( M,F ) = (0,0) (the primary modes), there exist secondary modes of instability 
for sufficiently large F. Instabilities of this type are shown in hgure[5] At M = 0, the first 
such secondary mode, which is odd, appears close to F = 2.3; the next secondary mode, 
which is even, appears close to F = 3.8. For larger values of F, the secondary modes 
are more unstable than the primary modes; this behaviour persists for M ^ 0, as shown 
in figure [5j>. Although an asymptotic description of the secondary modes is possible as 
F oo, the growth rates are relatively small, so these modes are not discussed further. 
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Figure 5. Instability of the rectangular jet at a — 1, and M = 0 (left) and F = 20 (right). 
Shown is c,' = Irn(c) computed numerically from (4.151 and (4.161. Even modes are shown as 
solid lines, and odd modes as dashed lines. The unique even and odd modes that smoothly 
connect to the modes at (A/, F) = (0,0) are shown as thicker lines. 


5. Smooth profiles: hyperbolic-tangent shear layer and Bickley jet 

As a model of more realistic velocity profiles, we consider in this section the instabil¬ 
ities of the hyperbolic-tangent shear layer and the Bickley jet. These may be regarded 
as smooth versions of the piecewise-linear profiles studied in (jdj Linear instability cal¬ 
culations involving these two profiles are well documented in a wide variety of contexts 


(e-g-, 

Lipps 

1962j Howard 1963 Michalke 

1964 

Sutherland & Peltier 

1992 Hughes & 

|Tobias|2001 

i and these provide a comparison and 

check on our results. We again restrict 


attention (for simplicity) to the case of a uniform background magnetic field, B(y) = 1, 
and with no underlying topography. 


5.1. Numerical method 


We seek a numerical solution of the eigenvalue equation (2.7), written as 


G 


S 2 


(K 2 y 

K 2 


G' - a 2 K 2 G = 0, 


(5.1) 


where, again, S 2 = [U — c) 2 — M 2 and K 2 = 1 — F 2 S 2 . Although the velocity profiles 
are defined over the entire real line, we solve (5.1) on y £ [—L,L\, using a shooting 


method, with matching imposed at y = 0, employing a generalised Newton method as 
the root-finding algorithm. Since the solutions decay exponentially as \y\ becomes large, 
namely 

G ~ exp (— aK±y) as \y\ —i oo, (5.2) 

where 

= l-F 2 ((U ± -c) 2 -M 2 ), U± = lim U(y), (5.3) 

y—>±oo 


K± 


we adopt expression (5.2) as the boundary condition to be implemented at y = ±L. Since 


our interest is in instabilities, singularities in the governing equation are avoided. To en¬ 
sure negligible influence of the finite domain, L is doubled until the computed eigenvalue 
changes by less than 0.5%. The routines are written in MATLAB, using odell3 as the 
integrator (an Adams-Bashforth method with adaptive grid). Although the boundary 
conditions are functions of c, changing at every iteration, we generally have no problems 
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with convergence provided that the initial guess is close to the true value. Solutions are 
initialised from (M, F) = (0, 0) at some fixed a using a known numerical result docu¬ 


mented in, for example, Drazin & Reid (1981). Runs at new parameter values are then 


initialised using an estimate for c from previously calculated values at nearby parameters. 
The Bickley jet is even about y = 0 and hence the parity result of §3.3| holds — i.e. the 
eigenfunctions are either even or odd. In this case we need integrate only up to y = 0, 
with the imposition of either G'(0) = 0 (even mode) or G(0) = 0 (odd mode). 


5.2. Hyperbolic-tangent shear layer 
In this subsection, we consider the basic state velocity defined by 

U(y) = tanhy. 


(5.4) 


From inequality (3.2) we know that the growth rates aci are bounded above by |I/ , | max /2 = 


1/2; furthermore, from the stability criteria (3.10) or (3.11), this profile is stable when 
M > 1. 

When M — F = 0, instability exists only when 0 < a < 1, with a neutral mode at 
a = 1 (e.g., Drazin & Reid 1981 §31.10). in this case, there is a single primary mode 


of instability, which may be classified as an inflection-point instability and attributed 
to interacting Rossby waves supported by the background shear (see, for example, the 
review by Carpenter et al. 2013). For both two-dimensional compressible hydrodynamics 
and shallow-water hydrodynamics, there is a secondary mode of instability, first found 


by Blumen et al. (1975). This has a smaller growth rate and a less pronounced spatial 
decay than the primary mode. Further, whilst the primary mode has c r = 0, there 
are two branches for the secondary mode with equal and opposite (non-zero) phase 
speeds, consistent with the parity results in §3.3[ which ensure that c = ±c r + i Ci for 
unstable modes. The secondary mode can be attributed to interacting gravity waves (e.g., 


Satomura|1981 Hayashi fe Young|1987 Takehiro fe Hayashi|1992 Balmforth| l999| ) and, 

indeed, can occur for linear shear flows, explicitly filtering out the possibility of Rossby 
waves due to a background vorticity gradient. This is consistent with the theorem of |Ripa| 
(1983), which states that for instability, either the associated potential vorticity profile 
possesses an inflection point or F > 1, both of which encourage instability through 
interacting waves. 

Both modes have been found in our SWMHD system. Figure [6] shows contours of Cj 
over ( M,F ) space at selected values of a for the profile ( |5.4[ ), distinguishing between 
the primary mode (with c r = 0 , shown as solid contours) and the secondary mode (with 
c r 7 ^ 0, shown as dashed contours). In figure^, at a = 0.7, and figure [ 6 j?, at a = 0.44 
(which is the most unstable mode when M = F = 0 ( Michalke||1964 [)), only the primary 
mode exists. In figure 0 ^ at a = 0 . 2 , and figure [ 6 ji, at a = 0 . 01 , both modes exist, 
although largely in different parts of (M, F) space. Note that figure 01 and figure [l] 
are remarkably similar, suggesting that long-wave instabilities for this velocity profile 
resemble those of the vortex sheet. This will be quantified via a long-wave asymptotic 
analysis in §5.5| in which we derive the more general result that long-wave instabilities 
of any shear-layer profile display the characteristics of vortex sheet instabilities. 

The growth rate aci is shown in figure [7] as a function of a. It can be seen that the 
secondary modes generally have weaker growth rates than the primary modes, consistent 


with the results of Blumen et al. (19751. As we shall see in f 5.5 the relation between the 


two types of unstable modes can be explored in some detail in the long wavelength limit. 
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Figure 6. Contours of a over ( M,F ) space at (a) a = 0.70, (6) a = 0.44, (c) a = 0.20, ( d ) 
a = 0.01. The primary modes (here classified by \c r \ < 10~' s ) are contoured as solid lines. The 
secondary modes are contoured as dashed lines. 



Figure 7. Growth rate versus a at selected parameter values for U(y) = tanhj/. The upper 
four curves are primary modes, the lower two are secondary modes. 


5.3. Instability mechanism: Counter-propagating Rossby Waves 

As mentioned above, inflection point instabilities can be attributed to interacting Rossby 
waves supported by the background shear. The constructive interference of a pair of 
Counter-propagating Rossby Waves (CRWs) has been proposed as the mechanism leading 
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Figure 8. Modified CRW mechanism in schematic form. Shown are two displaced material 
contours. The associated vorticity anomalies when M = F = 0 are shown by the closed solid 
curves; the effect of these on the other contours, leading to instability, is shown by the longer 
arrows. The closed dashed curves represent the additional vorticity anomalies when M and F 
are non-zero. The stabilising effect of these, which opposes the driving of the instability, is shown 
by the shorter arrows. 


to instability of shear flows in a variety of settings (e.g., Bretherton||1966 


Hoskins et al. 


1985 Baines & Mitsudera 1994 Heifetz et al. 2004| Harnik & Heifetz 2007). For the 
SWMHD system, it is therefore natural to enquire how this underlying mechanism is 
modified by magnetic and shallow-water effects. 

Let us first consider the case of M = 0, F = 0. The background vorticity profile sup¬ 
ports two Rossby waves, propagating in the negative (positive) ^-direction on the positive 
(negative) vorticity gradient in y > 0 (y < 0). Viewed individually, the Rossby waves 
are neutral and propagate against the mean flow. If, however, they become phase locked, 
they can interfere constructively, leading to mutual amplification and hence instability. 
This is shown schematically in figure [ 8 ] where the two Rossby waves are represented 
as perturbed vorticity contours (or equivalently perturbed material contours) for y > 0 
and y < 0. The resulting positive and negative vorticity anomalies are also shown. In 
this configuration, the transverse flow induced by each wave acts to amplify the existing 
transverse material displacement of the other. There is thus a mutual amplification and 
instability. This is also consistent with numerical solutions for instabilities of the flow 
(5.4) when M = F = 0; figure [ 9 J 2 , shows the vorticity perturbation of the most unstable 
mode (with a = 0.44), in agreement with figure [8] 

We now quantify how free-surface and magnetic effects modify this mechanism for 
SWMHD. We use the SWMHD vorticity equation, which is given by 

R '' Vw = -(V ■ u)w + M 2 B ■ S7j + M 2 (V • B)j, (5.5) 


Dt 


du> 

s m +u 


where u and j are the ^-components of the vorticity and electric current. Using equa¬ 


tions (2.1c) and (2.2), (5.5) can be written as 


Du; 

Dt 


= F 2 u-_ 


D h 
Dt 


+ F 2 huj'S/ ■ u + M 2 ( 1 - F 2 j)B ■ Vj - M 2 hjV • B. 


(5.6) 


On linearising about the basic state Uq 
of the form (2.3), and noting that v = 


= U(y)e x , Bo = e x , taking modal solutions 
(d/dt + Ud/dx)r], where y is the cross-stream 
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2ir/a 


Figure 9. Eigenfunctions of (a) vorticity and (6) height for the most unstable modes for 

U ( y) = tanh y at F = 0, M = 0. 


displacement, we obtain the the vorticity budget: 


u> = —r0 + F 2 hfl + M 2 


U-c 


(5.7) 


where SI = — U' is the basic state vorticity. The three contributions to w arise from the 
advection of the background vorticity and from shallow-water and magnetic effects. 

Inspection of figure [6] shows that the instability is most vigorous when M = F = 0; 
we therefore expect that the vorticity anomalies from the magnetic and shallow-water 


effects will be stabilising. The vorticity w and the decomposition (5.7) are shown in 


figure [10| for a mode at M = 0.25, F = 0.5. Even though this eigenfunction results from 
a calculation with non-zero F and M, the —rjQ 1 contribution has the same structure 
as that of figure The extra contributions from non-zero F and M are shown in 
figure [Tol u, d ; both of these terms are maximised at y = 0, where they are approximately 
in phase. Thus, at the simplest level, they lead to the vorticity anomalies shown by 
the dashed circles in figure [HJ The transverse flow induced by these vorticity anomalies 
counteracts the mutually amplifying transverse flow of the Rossby waves, and is thus 
stabilising. 

As a further verification of these ideas, we have also adopted a perturbative approach 


to the analysis of expression (5.7), approximating the shallow-water and magnetic con¬ 
tributions using the eigenfunction for M = F = 0. It can be seen that calculating F 2 hfl 
using h in figure [9j> is consistent with the full linear equations (figure [To] ^). To obtain an 
estimate of the magnetic contribution, it is necessary to calculate j using the governing 
equations ( |2.4[ ) with the velocity obtained when M = F = 0. This is slightly more in¬ 
volved than for h , but can be shown to provide a vorticity contribution consistent with 
figure \w\l (see Mak 2013). 


5.4. Bickley jet 

In this subsection we consider the basic state velocity defined by 

U{y) = sech 2 y, 


(5.8) 


again with B(y) = 1. From inequality (3.21, the growth rate aci is bounded above by 
I U r |max/2 = 2/(3\/3); furthermore, from stability criterion (3.11), this flow is stable when 


M ^ 1/2. When M = F = 0, even and odd modes are unstable only in the respective 
bandwidths 0 < a < 2 and 0 < a < 1 (e.g., Drazin & Reid||1981 §31.9). 

Figure [lT] shows contours of c* over (M, F) space for selected values of a. As for figure[6j 
the eigenvalues are calculated using a mode-tracking procedure, starting from the case of 
M = F = 0. The values of a correspond to: (i) the most unstable mode when M = F = 0 
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Figure 10. A generic vorticity budget breakdown for the case where neither F nor M is zero; 
shown here is F = 0.5 and M = 0.25, with (a) lo, ( b) —rjQ', (c) F 2 hfl, ( d ) M 2 j/(U — c). Notice 
that the vorticity contribution from the magnetic term is substantially larger than that from 
the shallow-water term. 


(panels a and 6); (ii) the mode with highest c* when M = F = 0 (panels c and d); (iii) a 
long-wave disturbance (panels e and /). In all these cases the magnetic field is stabilising. 
This effect will be quantified later via a long-wave asymptotic analysis. 

Figure[X2|shows the growth rate at selected parameter values. In general, the even mode 
is more unstable than the odd mode. There are isolated regions where the odd mode is 
more unstable, although these do not necessarily correspond to the regions predicted by 
the stability analysis for the rectangular jet, described in §4. 

A natural question to ask, prompted by the findings for the shear layer as well as 
the results for the rectangular jet in §4.2.5, is whether there are additional modes of 
instability to those shown in figure [l2] We have performed a scan over (M, F ) space at 
various values of a , with randomly generated initial guesses for c within the smallest 
rectangle containing the semi-circle (3.9). A substantial number of computations (20 
different initialisations at over 200 different parameter values) were carried out, solving 
the governing eigenvalue equation with no parity imposed. Other solution branches were 
found for large F when a is sufficiently large, analogous to the secondary modes for the 
rectangular jet in §4.2.5. Sample plots of c, and cross-sections of the eigenfunctions for 
the even modes are shown in figure |13[ the results for the odd mode are similar. As in 
§4.2.5, further secondary modes appear as F increases (figure 13 z) and the associated 


growth rates are small (figures |13| z, b). However, in contrast to the results shown in 
figure [5] for the rectangular jet, we did not find branch crossings in c,; for the Bickley jet; 
this may occur at higher values of F. but computations are demanding since c* becomes 
increasingly small. 

There are two distinctive features of the spatial structure of the eigenfunctions, as 
shown in figures [13} :, d, e. First, there is a core region consisting of approximately quan¬ 
tised oscillations: the primary mode has one oscillation, the n th secondary mode has 
n + 1 oscillations. Second, the boundary of the core region is approximately located 
where U(y) = c r . Here the eigenfunctions are highly oscillatory and of larger amplitude, 
although they are well-resolved in figures [Ki}:, d, e. 
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Figure 11. Contours of d over (M, F) space at selected a for the even (left column) and odd 
(right column) modes of the Bickley jet, with (a) a = 0. 9, (b) a = 0.52, (c) a = 0.3, ( d ) a = 0.23, 
(e, f) a = 0.05. The cutoff from the asymptotic result (5.241 is also plotted in panel (e). 



Figure 12. Growth rate of the primary instability of the Bickley jet. 


5.5. Long-wave asymptotics for unbounded smooth profiles 
Further understanding of shear flow instabilities in SWMHD can be obtained by gener¬ 
alising the long-wave asymptotic procedure of D razin fe Howard] (1962), who considered 
two-dimensional hydrodynamics. The idea is that, for long-wave disturbances, the lead¬ 
ing order behaviour of the instability is determined by U(y) as \y\ —> oo, with higher 












































Figure 13. Even mode of the Bickley jet at a = 1. (a) d as a function of F at M = 0, a = 1. 
(6) d as a function of M at F = 20. Crosses show the primary mode, circles and triangles the 
first two secondary modes, (c, d, e) Cross-sections of the vorticity perturbation for M = 0.05, 
F = 20 for the primary and first two secondary modes. The locations where U = c r are labelled 
on the axes. The spikes reach an amplitude of approximately 0.5; the display range has been 
reduced in order that the oscillations within the core are clearly visible. 


order corrections determined by the flow at finite y. In this subsection we extend the 
formalism of Drazin & Howard (1962) to SWMHD, but with a uniform magnetic field 
and no topography. 


We consider the governing equation (5.1), written as 
Z 2 {G" - a 2 K 2 G) + ( Z 2 )'G' = 0, 


Z 2 = 


K 2 ' 


(5.9) 


We assume that U± = t/(±oo) are well defined and that U' (and so (Z 2 )') decays 
sufficiently rapidly as \y\ —»• oo. Then, on choosing an appropriate frame of reference and 
suitable normalisation for the basic flow, any velocity profile may be designated as either 
a shear layer if U± = ±1, or as a jet if U± = 0. 

Adopting the same notation as Drazin & Howard (1962), we consider solutions to (5.9) 
of the form 


G(y) = 


G+(y ) = x(y) exp (-aK+y), y > 0, 
G-(y) = 9{y) exp (+aK_y), y < 0, 


(5.10) 


where K± = 1 — F 2 S’± = 1 — F 2 (( U± — c) 2 — M 2 ). The perturbations must decay 
\y\ —► oo; hence Re(/C±) > 0. We consider expansions of the form 


as 


x(y) = ^2(+a) n Xn(y), 9{y) = ^2(-a) n d n (y), 

n =0 n —0 


(5.11) 
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with Xo,8o —> constant 0) and Xn,0 n —> 0 as \y\ —> oo (n G 1). It turns out to be 
most convenient to fix xo(°°) = $o(— oo) = 1 , and then to accommodate the necessary 
degree of freedom in the matching conditions for G at y = 0, namely G+(0) = FG_(0) 
and G' + (0) = TG(_( 0 ) for some constant F. Consistency thus implies 

G+( 0 )G'_( 0 ) = G_(0)G / + (0). (5.12) 


Without loss of generality, we shall focus on the equations for those for 0 follow in a 
similar fashion. On substituting (5.10) (with ( |5.11| )) into (5.9), equating the coefficients 
at each order of a gives 

0=(Z 2 x' o Y, (5.13a) 

0 = (zVi)' - K+ (2ZY 0 + ( Z 2 )' X0 ), (5.136) 

o = {Z 2 Xn+ 2 )' - K + (2ZY n+ i + ( z2 YXn+ 1 ) + Z 2 (I<1 - K 2 ) Xn , n > 0. (5.13c) 

Equation ( |5.13a| ) integrates to Z 2 x'o = G, with the conditions at infinity then giving 
C = 0. Thus xo = constant = 1 through our choice of xo(oo)- Integration of equations 
(5.13f>, c) then gives, after some algebra, 


Xi = 


X 2 

1 —+ 

z 2 


dyi, 


» " II 0 /> 2 ~ sl) dw+ ** £ ( d - ¥) d ») d *" 


(5.14) 


The matching condition (5.12) then leads to the result 

S 2 


0 = 


4 

K, 


KG 


( /* OO nO 

/ ( S 2 ~S 2 + )dy+ (S 2 -S 2 _)dy 

JO J- oo 


(5.15) 


4 


K+K- 


1 - 


S 2 


d y - 


K+KG 


!- ^ ) dy) +0(a 2 ). 


equation (5.15) then 


On expressing the eigenvalue c as c = + ac^ + a 2 c^ + 

determines the successive cY 

Although we have focused on the case of a uniform magnetic field, it is possible to 
include a non-uniform field, subject to imposing conditions analogous to those for U{y). 


With underlying topography, other assumptions are required (Codings & Grimshaw 


1980). 


5.5.1. Shear layers 


For a shear layer, U± = ±1, and the leading order term of expression (5.151 gives 
(1 - C (°>) 2 - M 2 


(1 + C (°>) 2 - M 2 
yjl - F 2 ((1 - c(°)) 2 - M 2 ) \Jl~F 2 ((1 + c(°)) 2 -M 2 ) 

This is exactly the eigenvalue equation of the vortex sheet 
layer, c —> Cv, as defined in (4.7), as a 


= 0 . 


(5.16) 


4.5[ ); hence, for any shear 
0. This is not surprising; sufficiently long waves 
see the shear layer as a discontinuity in the flow. 

When F is sufficiently large, c v is real (see §4.1). In this case, following Blumen et al. 


(1975), we calculate to seek a secondary mode of instability. For the particular case 
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U(y) = tanh y. and after considerable algebra, we obtain 

y/F 2 ((1 - c v ) 2 - M 2 ) - 1 (1 + c v ) 2 - M 2 
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4c, 2 , 


1 + f ( log 


(1 


y/l + 4E 2 + 4 F 4 M 2 
c v ) 2 - M 2 


(1 - c v ) 2 - M 2 


- - 6 , 


(5.17) 


1 - M 2 - i 
4M 


log 


(1 + M) 2 - eg 
(1 - M) 2 - c 2 


+ < 5 i — 82 


where 


61 = 


0 , 


|c v -M | < 1, 
otherwise, 


= 


0 , 


c v + .M ^ 1, 
otherwise. 


It may be shown that for M —>• 0, the expression inside the square brackets is equal 
to 4M(1 — c v ) -2 + 0(M 2 ), and that equation (5.171 is equivalent to equation (21) in 


Blumen et al. (1975]). Equation (5.171 does indeed describe instability beyond the vortex 
sheet cutoff. By expanding c v up to powers of it may be shown from (5.17) that 
In^c*- 1 )) \, 0 as F —> 00 , so there is no cutoff at finite F. 

The analysis leading to expression (5.17) is valid only when c v is not small, i.e. for F 2 
not close to 2(1 — M 2 ) -1 or M not close to 1. When 2 — F 2 ( 1 — M 2 )~ l is small, 0(a 2 / 3 ) 
to be precise, we have c v ~ ac^ ~ a 1 / 3 . Rescaling and choosing the appropriate branch 
so that Re(v /rTT ) > 0, gives c 11 ^ as the solution of the cubic equation 

6 + 2M 2 / m\ 3 / 2 ni . ( 1 


(1 - M 2 ) 3 




1 - M 2 


- F‘ 


c*- 1 ! — ia 


1 - M 2 


1 , 1 + M 

2 M ° g 1 - M 


= 0, (5.18) 


for a 2 / 3 ~ 2(1 — M 2 ) 1 — F 2 —> 0. When M = 0, equation (5.18) reduces to equation (23) 
Blumen et al. (1975]). There are two admissible roots with positive imaginary parts. 


There is a transition to non-zero real parts when 


F 2 = 


cusp 1 - M 2 


- 3 


a 2 6 + 2M 2 
T (1 - M 2 ) 3 V1 ~ M 2 


1 1 + M' 

2M ° S 1 - M. 


1/3 


(5.19) 


This expression reduces to E 2 usp = 2 — 3(6a 2 ) 1 / 3 when M 


0 , as given by 


(1975). 


Blumen et al. 


The asymptotic results, together with the numerical computations of §5.2, are pre¬ 
sented in figure |14| The agreement between the two, including the location of the cusp 
given by equation (5.19), is excellent. 


5.5.2. Jets 


For a jet, U± =0 and (5.15) simplifies to 


0 = 


2Sg 

Kn 


/ (5--*)*,-3$/ ( 


1-f) d y 


0(c 


(5.20) 


where Sq = (0 — c) 2 — M 2 . Here, for a fixed value of F, we need to consider different 
regimes for M. 

For F 2 = 0(1), if M 2 = 0(1) then c*- 0 -* and c^ are real. To find an instability we 
need to consider the regime M 2 ~ a, which implies c^ 0 ) = 0. At the next order, we 
choose to balance the first two terms on the right hand side of ( |5.20[ ), assuming that the 
second integral is negligible; this is confirmed by the analysis of Appendix [A] Defining 
E = U 2 /2 d y, where U is assumed to decay sufficiently rapidly that E is finite, we 
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Figure 14. a and c r versus F at a = 0.01 for (a) M = 0, (6) M = 0.25, (c) M = 0.5. The 
crosses are the computed results for U(y) = tanhy, the solid lines are the asymptotic result 


c v + qc^ 1 -* with f' 1 - 1 given by (5.171, and th e dot- dashed lines are the inner expansion given by 
the relevant solution to the cubic equation (5.181. 


obtain 

c ~ iy/aE - M 2 as M 2 ~ a —>• 0, F 2 = 0(1). (5.21) 

The corresponding result for compressible hydrodynamics was derived by Gill & Drazin 
(1965), and for incompressible MHD by Gedzelman (1973). 

For large F, the regime of interest is M 2 ~ a, F 2 ~ a~ L . Considering the same balance 
as above gives 

c ~ i - M 2 + v '° 2£2 + 1/2 as 0. 

(5.22) 


This result reduces to (5.21) in the limit of small F. Equations (5.21) and (5.22) are 


the extensions to smooth velocity profiles of equations (4.24) and (4.26), which are the 


dispersion relations for the even mode of the top-hat jet. It is interesting to note that 
the odd modes are not recovered by this analysis. 


For the Bickley jet, 2 E = (sech 2 y) 2 d y = 4/3 and (5.21) and (5.22) become 


c ~ i\l-a — M 2 


as 


and 


# 2 ^" 2 GV» 2 + vr 


M 2 


1/2 


0, F 2 = 0(1), 


(5.23) 


as 


F 


-2 


M 2 


a -» 0. (5.24) 


In figure 15 the growth rates given by ( |5.23| ) and (5.24) are plotted against those deter¬ 
mined numerically from equation (E3; there is good agreement at small a. The cutoff 
implied by (5.24| is also shown in figure [TT^ . 

The above asymptotic procedure does not yield the odd modes for smooth velocity 
profiles; for these modes, c r remains 0(1) as a —> 0, as can be verified numerically. 
The difficulty can be traced back to equation (5.9) since, when M 2 ~ a, Z 2 becomes 
small when U « c r . The standard asymptotic procedure leading to equations (5.13) 
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a 


Figure 15. Comparison of the computed growth rates for the even modes (crosses) and the 
predicted growth rates from the asymptotic resul ts for U(y) = sech 2 y, for M = 0.1 and (a) 
F = 0, (b) F = 1, (c) F = 5, ( d ) F = 10. a x Im( 5.23j ) is given by the dot-dashed line (cutoff 
plotted as vertical dotted line) and a x Im(5.24| is given by the solid line (cutoff plotted as 
vertical dashed line). 


breaks down. Instead it is found that Ci ~ a 2 / 3 , 
hydrodynamic case by Drazin & Howard (1962). 


a result that has been derived for the 


6. Conclusions and discussions 


The SWMHD equations, introduced by Gilman (2000), are a useful model for studying 


MHD in thin stratified fluid layers. We have investigated the linear instability of paral¬ 
lel shear flows with an aligned magnetic field in planar geometry with no background 
rotation. The instability of hydrodynamic shear flows is a classical problem, with well 
known results such as Rayleigh’s inflection point criterion, Howard’s semi-circle theorem 
and Hpiland’s growth rate bound, supplemented by extensive asymptotic and numerical 
results for various idealised flows. Motivated by geophysical and astrophysical consider¬ 
ations, previous authors have extended separately the analysis to include the influence 
of shallow-water dynamics and a magnetic field. Here, for the first time, we have applied 
this classical approach to the SWMHD system, complementing previous work on the 
instability of specific flow configurations in spherical geometry ( Gilman &; Dikpati||2002 
Dikpati et aL|[2003 ). 


We first considered the stability of arbitrary flow and field profiles, leading to the 
growth rate bound (3.2) and the semi-circle theorems (3.7) and (3.9), which, in turn, 


imply the stability criteria (3.10) (3.12). More detailed investigations focused on shear 

layers and jets, which are known from hydrodynamical problems to be canonical flows 


(Drazin & Howard 1962). For each of these we considered both discontinuous profiles, for 


which dispersion relations can be obtained analytically, and analogous smooth profiles, for 
which numerical solutions were obtained. We also took the background magnetic field 
to be uniform; this simple choice already leads to interesting stability characteristics. 
The imposed field strength is measured by a parameter M, whilst the stratification is 
measured by the Froude number F. 
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For the shear layer we considered the vortex sheet and the hyperbolic tangent velocity 
profile. For both profiles, a key finding is that, although increasing M or F in the absence 
of the other is seen to be stabilising, their combined effects can offset each other, resulting 
in a tongue of instability for arbitrarily large F and M approaching unity. However, the 
strongest instabilities are found at smaller F and M. For the smooth shear layer, these 
were interpreted through the counter-propagating Rossby waves (CRWs) mechanism, 


which has been used to understand hydrodynamic shear instabilities (e.g., Carpenter 
et al. 2013). Here we show that the vorticity anomalies associated with shallow-water 


effects and magnetic tension oppose those for the underlying CRW mechanism. As far as 
we can tell, this is the first time such a dynamic argument has been invoked to explain 
the stabilisation of shear flow instabilities by a magnetic field. For the smooth shear layer 
alone, there is a weak secondary instability when F > 1, analogous to those found in 
other hydrodynamic systems ( Blumen et al. |[l975 Balmforth||1999 ). Further insight was 
obtained by performing a long wave asymptotic analysis for an arbitrary shear profile, 


extending the hydrodynamic treatment of Drazin & Howard (1962) to SWMHD. There 


is a primary mode of instability, yielding a dispersion relation identical to that of the 
vortex sheet, thus showing the widespread applicability of the vortex sheet model. For 
sufficiently large F, the long wave analysis recovers the weaker secondary mode. 

For jets, we considered the top-hat velocity profile and, as its smooth counterpart, 
the Bickley jet. For both profiles, the modes of instability are either even or odd about 
the jet axis. At large a , dynamics on the flanks of the top-hat become uncoupled, and 
hence the instabilities become those of the vortex sheets on the flanks. However, this 
behaviour does not carry over to the Bickley jet, since the even and odd modes are 
stabilised when a > 2 and a > 1, respectively. For small a, the even mode is only 
unstable when M < M c = 0(a 1//2 ). This result applies for both the top-hat profile and 
arbitrary smooth jets, showing the widespread applicability of the top-hat model for the 
even mode at small a. For small a, the odd mode is also only unstable for sufficiently 
small M. The growth rate Cj = 0(a 1 / 2 ) for the top-hat profile, but Ci = 0(a 2 / 3 ) for 
the Bickley jet, so the analogy between the top-hat and smooth velocity profiles is lost. 
Although we can be precise about the behaviour at large and small a , in general either 
the even or odd mode may be more unstable. Finally, in addition to these two primary 
modes of instability, there exist multiple secondary modes of instability at large F, with 
relatively small growth rates. 

Given the motivation for this study, it is important to consider the implications of our 
results for the solar tacliocline. Adopting the tachocline radius (« 0 -7Rq) as a character¬ 
istic length scale L 0 and taking Uq = 500ms -1 as a characteristic velocity, as discussed 
in §2, the e-folding time for an instability is given by t = ( aCi)~ 1 (L 0 /Uo ) « 10 6 (aCi) -1 s. 
With F ss 0.1, aCi ss 0.1 over the range of M anticipated for the tachocline, for both shear 
layers and jets (see figures [?| and [l2|) , leading to large-scale instabilities with t = O(10 7 )s, 
a relevant timescale for solar dynamics. 

It is hard to draw direct comparisons with previous results on SWMHD instabilities 
( Gilman fe Dikpati|2002 Dikpati et al. |20f)3 Gilman &; Cally|2007 ). These studies, which 
were focused on the tachocline, were in spherical geometry with background rotation. Fur¬ 
thermore, the main emphasis of the work was on the destabilisation of hydrodynamically 
stable velocity profiles by non-uniform toroidal fields. 

Having investigated the quite complicated behaviour that results from our idealised 
model, there are, within the present geometry, several natural extensions of our study. The 
inclusion of rotation could be considered; even in the hydrodynamic case this modifies 
the nature of the instabilities. Non-uniform magnetic fields open up the possibility of 


destabilising hydrodynamically stable velocity profiles, as in Gilman & Dikpati (2002) 
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and Dikpati et al. (2003). The nonlinear evolution of the instability and the associated 


changes in the mean flow are clearly of interest; some preliminary results are given in 


Mak (2013). 


This work was supported by the STFC doctoral training grant ST/F006934/1. JM 
thanks Eyal Heifetz for helpful discussions. 


Appendix A. Consistency checks for the long-wavelength jet analysis 

The aim of this appendix is to show that 


c 2 - M 2 


p+oo 


1 — F 2 (c 2 — M 2 ) J_ 0 


1 - 


c 2 - M 2 
(U - c) 2 - M 2 


d y 


(Al) 


is O (o;log( a)), and hence that the a symptotic analysis of § 5.5.2 is consistent. 

Following Drazin & Howard (1962), we assume that \U\ ^ Ae _ °l y l, which is satisfied for 
the Bickley jet. They make the additional assumption that c is ‘almost pure imaginary’; 
here we adopt the modified assumption that 




which is supported by both our numerical and asymptotic results. 

d y. 


(A 2) 


(A3) 


(A 4) 


Consider first the case of F = 0 for which (All is given by 

' ^ U 2 ) ./ I .( 1 (C-c) 2 - \ 1 \ 

For y > 0, we split the range of integration into (0, A) and (A, oo), where 

/ A\ 1/a 

A = log (—J = O (log(a)), 

since c 2 = 0(a). Then 

(c 2 - M 2 ) J (■ ■ •) dy < |c 2 - M 2 \£ (l + 

< |c 2 -M 2 |A(1+A) = 0(a log a), (A5) 

where we have used the usual integral inequalities, the inequality \(U — c) 2 — M 2 | A c 2 , 
the assumption that M 2 = O(a), and the derived result (5.21) that c = (^(a 1 / 2 ). 
Similarly, 


i cr + M 2 \ , 

2 , dy 

^2 i\/r2 


/ oo roc / 

{■ ■ ■) dy = (c 2 — M 2 ) [ 


< 


— M 2 \ 


2 , , , U 2 - 2Uc A 

a {(U- c ) 2 -FPJ dy 

(A 2 e~ 2aX : 2A|c|e" aA 


/ (U 2 + 2\U\\c\)dy^N ^ 


+ 


a 


= O(a). 


The dominant contribution is from (A5), and hence I = O(aloga). 


When F is non-zero, the only difference is in the pre-factor to the integral (Al). 


However, with c from either (5.23) or (5.24) the pre-factor is 0(1) and we may conclude 
that I = 0(a log a). 
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